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Abstract 

We discuss a phase space description of the photon number distri- 
bution of non classical states which is based on Husimi's Q{a) function 
and does not rely in the WKB approximation. We illustrate this ap- 
proach using the examples of displaced number states and two photon 
coherent states and show it to provide an efficient method for comput- 
ing and interpreting the photon number distribution . This result is 
interesting in particular for the two photon coherent states which, for 
high squeezing, have the probabilities of even and odd photon numbers 
oscillating independently. 
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1 Introduction 



The oscillations in the photon number distribution are an interesting feature 
common to various kind of light states, which may be taken as a signal of 
non classical behavior of those states. They were first computed for squeezed 
states P3 0] and interpreted as an interference effect in phase space Pfl El 
El 13]. They also have been the subject of experimental investigations [S] in 
connection with the properties of Wigner's distribution. Since in the actual 
experimental situation many factors like the detector properties may modify 
the conclusions of the theoretical analysis it is of interest to understand better 
the physics involved and to develop our intuition of how the oscillations may 
appear and eventually disappear. 

For our analysis we consider a single mode of the electromagnetic field. 
For this case, the picture of quantum interference of states in phase space was 
developed as a generalization of the Bohr-Sommerfeld description of quantum 
states as finite areas in phase space [H 121 El El Ej • The total area of each state 
is fixed by the normalization requirement {'4'\'4') = 1- The area associated 
with a Fock or number state |m > is a circular ring centered at the origin 
with interior radius r~ = v^2m and exterior radius r+ = ^/2/m+^. The areas 
associated to coherent states and squeezed states are obtained by displacing, 
or squeezing and displacing the one associated to the vacuum state. The 
inner product < > between two quantum states is then related with the 
overlapping area of the states in the phase space. The intersecting region 
may have in some interesting cases more than one component. Since the 
probability amplitude < ml?/' > is a complex number and the overlapping 
areas are real, it results natural to associate to each of these components a 
phase in order to reproduce the quantum mechanical results. The probability 
amplitude acquires then the following structure, 

< Xl^ >= E /4^e^^'^^^^ (1) 

i 

where A^^^{iIj) and (f>':^^{ip) are the i-th component of the overlapping area 
and their assigned phase respectively. 

The presence of the areas (■?/') in this expression is very natural and 
physically appealing although the fact that what appears is the square root 
of the areas does not allow any direct geometrical method to derive Eq.((T]). 
On the other hand the values that have to be chosen for the phases 
are not evident from the geometry of the pase space. Dowling [et a/.)0| 
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work out a quite general equivalent of Eq. ((T]) for the probability amplitude 
of the eigenstates of two different hamiltonian operators using the WKB 
approximation. They obtained an explicit expression of the phases where 
it is possible to recognize a geometrical content. Oscillations in the photon 
statistics for displaced states, two photon squeezed states and for squeezed 
number states may also be studied with this methodology but this approach 
is limited by the validity of the WKB approximation [Bl Ej • 

Some time after, Milburn jS] show that the interference effects in phase 
space may also be understood by considering the properties of Husimi's func- 
tion ^ 

Q{a) = -| < a\ip > p , 

TT 

for the state under study. This author notes that the over-completeness of 
the coherent states allows to rewrite the probability amplitude {mlip) as a 
phase space integral, 

< >= - J (fa < m\a >< a\ip > . (2) 

The functions < m\a > and < a\ip > could be interpreted as the phase space 
probability amplitudes for the states \m > and {ip > respectively. Each of 
these functions is proportional to the function Q{a) of the corresponding 
state. Then, 

<m\a>= v/7rQ„(a)e*'^'"(") (3) 

< a\tlj >= ^7rQ^(a)e'^*(") (4) 

One may then approximate the integral in (0), by restricting the domain of 
integration to the phase space region where the product of the two probability 
amplitudes is appreciably different of zero and by identifying the regions were 
the Q function concentrates with the phase space regions assigned to the 
states, this approximation is equivalent to consider the integration domain 
as the overlapping areas between the different states. Following this line of 
thought in this paper we discuss in detail the photon number distributions 
for the displaced number states and the two photon coherent states. We 
show that this approach allows to identify the areas and phases for the phase 
space description without rendering in the WKB approximation. This results 
is interesting in particular for the two photon coherent states which has, for 
high squeezing, the probabilities of even and odd photon numbers oscillating 
independently. 
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2 Photon statistics for displaced number states 

Let us consider first the example of tlie displaced number states. They are 
defined through the action of the displacement operator D{j3) = exp{Pa^ — 
13* a) on the Fock states |n), 

\n,[3) = Dmn) ■ (5) 

For simplicity we work with real j3. 

We compute the photon distribution using the method of generating func- 
tions. Consider a coherent state \a) = D{a)\0) with real a. We have, 

n=0 

Then, 
with 

1 1 2 

{m\a + /3) = ^=e~^^ {a + /?)'". (8) 
V ml 

As in other cases the photon number distribution PmniP) = P)]"^ is 

oscillating, (see Figure (jHl) below) 

To develop the phase space description consider the phase space ampli- 
tudes, 

{m\a) = ^e-l^l/'a™ 
vm! 



and 



|(m|a)| e''^'"(") (9) 



1 I3a*-al3* \a-fJ\^ 

(a|n,/3) = ^e^-^e-^{a* - (3T 
vn! 

= |(a|n,/3)| e-*^"-''^"). (10) 

The phases here are better expressed in terms of the real an imaginary parts 
of a = X + iy and take the form 

0m(a) = marctan ^— j (11) 



and 




Figure 1: | < m|a > | for a Fock state with m = 3 

In FigureilH) we show the phase space probabihty amphtude module 
|(m|Q;)| for a number state with m = 3. This function reminds the Bohr- 
Sommerfeld ring associated with the corresponding number state except for 
the fact that here the mean radius is |a| = -y/m. This suggests to associate 
the number state |m) the circular ring centered at origin and located between 
the radii \/ m— 1/2 and + 1/2. Note also that the phase space ampli- 
tude |(a|n, /?)| for the displaced number state |n, /?) is obtained by displacing 
the amplitude |(a|n)|. 

In Figure Q we show the product | < m|a >< > | for m = 

100, n = 3 and (3 = 10.1. Observe that there are two overlapping regions 
where both terms are appreciably not vanishing. But note also that the 
intersection areas with this prescription are not in the same position were 
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Figure 2: | < m\a >< a\n,l3 > \ for a Fock state with m = 100 and a 
displaced state with n = 3 and (3 = 10.1 



they would appear if working with the Bohr-Sommerfeld bands. Now they 
may be localized in the intersection points of a circumference of radius ^/m 
centered at origin and a second circumference of radius y/n displaced by (3, as 
shown in Figure©. The intersection points are then given by a;+ = Xq + iyo 
and = Xq — iyo, with 

f3^ + m-n 
^0 = 2/3 

and 



y^=^m-xl (14) 

Now we approximate the phase space integral of Equation ^ for the 
photon number probability amplitude {m\n, (5) by taking instead of Q and 
|TO|l uniform contributions on the intersection of the rings described above 
with the angles (|TT|l and (|T^ evaluated at the points given by (|T!^ and (fT^ . 
The probability amplitude has then the structure, 



(15) 



Figure 3: Intersection points in phase space 




Figure 4: Overlapping areas between the Fock state \m > and the displaced 
number state \n, j3 > 



The computation of this areas is analogous to the one presented in j^l E] 
for the Bohr-Sommerfeld strips and is given by, 

Amn = \ {a{^m + 1/2, + 1/2) - a{^m - 1/2, + 1/2) 

-a{^m + 1/2, - 1/2) + a{^m - 1/2, - 1/2)) (17) 



where a(r, i?) = + i?^7 — /JYi is the internal area between the two circular 
paths with radii r and R as shown in Figure (jSj) , with (Xi = ^ "'"2/3"^ ? ^1 = 






B 



Figure 5: Internal areas used to define A^^ 

Figure (0) show the exact and approximate results of Pmn for a displaced 
state with = 3 and /3 = 10.1 showing a good qualitative and quantitative 
behavior for the values for which the assigned rings actually overlap. 

The comparison of the result in Eq. (fTBj) which has the direct geometrical 
interpretation 

V'S = a(v^, Vn) - nir (18) 

with the one obtained by Dowling et al using the WKB approach is shown 
in Figure ((Tj). 
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Figure 6: Photon statistics Pmn for a displaced number state \n, P > with 
n = 3 and (3 = 10.1. Exact computation: crosses. Approximationxircles. 



The WKB phase is given by gj E] 



WKB 



(m + 1/2) arcsin 



Xc['m, n) 



+ (n + 1/2) arcsin 



V2m+ 1. 

a/2/3 / -TT in — m)n n 

V2m + 1 - Xc[m,nY ^ h - 



Xc{m,nj 



V2n + 1 



(19) 



where 



m — n \/2(3 
Xc[m,n) = — + 



■ 2 (2°' 

The phases differ sensibly for zero and high photon number where the WKB 
approximation fails to reproduce the distribution [7]. 



3 Two photon coherent states 

As a second illustration let us consider the Two photon coherent states. 
These are obtained by the application of the squeezing operator S{r) = 
exp — |('^^)^} on the coherent state That is, 

|/3,r) = 5(r)|/5) (21) 
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Figure 7: Comparison between the phases. Crosses: {if) + mi) / ti) . Circles: 



Again we take (3 and ^ real. The photon statistics is given in this case by 

unmn, 

(tanh(r))' 



P 



\{n\P,r)Y 




exp {/5^(tanh(r) - 1)} 



2 cosh(r) sinh(r) 



(22) 



with Hn{x) the order n Hermite polynomial. 

For r << 1 the photon statistics resembles that of a coherent state (al- 
though in fact the statistics is sub-poissonian) . As the squeezing increases, 
they appear the oscillations in the distribution characteristic of the phase 
space interference. But then for higher values of r ( see Figure (jSI) ) the 
distribution develops different oscillating behaviors for odd and even photon 
numbers an effect which, as we show below, may be understood in terms of 
phase space interference. Following the same line as in the last section take 
Equations Q and (fTT|) and consider. 



(a|/3,r) 



1 



a/ sech(r) exp <^ — -(|ap + /3^) + a*/3sech(r) 



-i((a*)2-/32)tanh(r) 
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Figure 8: Photon statistics for a two photon coherent state with /3 = 5.1 and 
r = 3. 

= |(a|^,r)| e''^''-(") (23) 

which define 

0/3,r(«) = — ?//3sech(r) +x?/tanh(r) . (24) 

It is not difficult to show that the phase space probabihty amphtude 
|(q;|/3, r)| of the two photon coherent state concentrate on an elhpse centered 
in /3exp{— r}, with one semi-axis of magnitude e~'^ oriented through the X 
axis and the other semi-axis of magnitude e*" oriented through the Y axis. 
We represent the state by the internal region of the ellipse, 

(x-/3e-'-)V'^ + yV2'- = 1. (25) 

of area vr. For high values of r the ellipse approximates to a vertical line. 
In this limit the intersection points with the circumference of radius ^/rn 
(which we are using to represent the number states) are a+ = X2 + iY2 and 
a_ = X2 — iY2, with 

X2 = /?exp{-2r} 

Y2 = ^Jm-Xl . (26) 
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Figure 9: Intersecting points which the phase for < m|/3, r > 



As in the previous case we may write the probability amplitude (m|/?, r) 

as: 

(^1/5, r) = ^e^'^^- + ^^e-'<r (27) 

TT TT 

where the phases evaluated at the intersection points are given by 

. (m) 



^13 J = +0/3,r(a+) 
f Y2 1 

= m arctan < — > — Y2Psech.{r) + X2Y2 tanh(r) (28) 

1^2 J 

and Am is one the shadowed areas in Figure ()10|) . It has the value 

Am = AyAx (29) 

with 

Ay = v/m + 1/2 -/52g-2r _ _ x/2 - /?2e-2r 

Ax = 2e-Vl - (m-/52e-2'^))e-2'^ 

Figures (lllj) and ()12|) show the comparison between the exact and ap- 
proximate results for the photon statistics of a two photon coherent state. 

To get a geometrical interpretation of the two branches in the photon 
distribution consider again the phase evaluated at (j2Hl)- It may be rewritten 
in the form, 

^ = m0 - X2Y2 (^j^ - ^^^^] =m<p- X2Y2 (30) 
\cosh(r) cosh(r) / 
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Figure 10: Overlapping for < m\(3, r > 



From this expression follows that the phase V'/3,r is represented by the shad- 
owed area in Figure Q- For high values of r this phase may be written 
as, 

^lJf,,r = m7i/2 - 2X2Y2 (31) 
In this case the probability is given by, 

P^ = — cos\^Pp,r) = 4^ [1 + (-1)™ cos{AX,Y2)] (32) 

which explains the differences in the probabilities of odd and even photon 
numbers. 



4 Conclusion 

In the approach to quantum phase space interference presented in this paper 
we show a general setup, based in the properties of Husimi's Q{a) function, 
for the representation of quantum states as regions in phase space with a 
prescription for assigning phases to the different overlapping regions with 
a has geometrical input. We illustrate the method discussing the photon 
number distribution of the displaced number states and of the two photon 
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coherent states. For the latter at high squeezing the distribution we show 
that the distribution develops different oscillating behaviors for odd and even 
photon numbers. This effect (also displayed by squeezed states) is understood 
in terms of phase space interference and maybe, in principle, suitable for 
experimental verification. 
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Figure 11: Pm with m even for a two photon coherent state with /3 = 
and r = 3. Circles: Exact resuh. Crosses: Approximate result. 
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Figure 12: Pm with m odd for a two photon coherent state with /5 = 5.1 
r = 3. Circles: Exact result. Crosses: Approximate result. 
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